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Abstract 



This paper suggests a method to create dark soliton-like waves in a 
trapped atomic Bose-Einstein condensate (BEC). With the phase im- 
printing method, a soliton-like disturbance is produced along the x-axis. 
The method is used in two cases (with and without rotation of a trap- 
q ping potential). The disturbance in a two-dimensional (2D) system is 

£J , unstable. The soliton-like perturbation is found to decay into more sta- 

ble vortex pairs together with a burst of sound wave. The simulations 
*y-s \ are performed by solving the 2D Gross-Pitaevskii equation (2D GPE). 

^. ■ The successive dynamics of the wave function are studied by monitoring 

C*~) \ the evolution of density and phase profile. 

On 
00 

^ ! 1 Introduction 

r-> 

:rt: ■ Weakly interacting atomic gases, like BECs, give rise to study the non- 

linear properties of matter. Soliton solutions of the non- linear Schrodinger 
equation (NLSE) first were deduced by Zakharov and Shabat pp. Dynamics 
and stability of dark solitons were described theoretically [2J together with a 
method of their creation [3] . The effectively repulsive atomic interactions and 
the scattering properties between the atoms of the condensate support the 
appearance of solitons. 

If the scattering coefficient, g is positive (g > 0), the effective nonlinearity 
is repulsive and describes a dark soliton. Here, g = 4irh 2 a/m, a is the s-wave 
scattering length and m is the mass of an atom. 



Dark soliton in matter waves is characterized by a local density mini- 
mum and a sharp phase gradient of the wave function at the position of the 
minimum. Dark solitary wave is not a topological defect like vortices. 

In a dissipative environment, it can disappear anywhere in the conden- 
sate. Its quantitative study with a dissipative dynamics at finite temperature 
was well described [lj. At finite temperature, the dynamics of a vortex within 
the Zaremba-Nikuni-Griffin formalism was studied in atomic Bose-condensed 
gas and clarifid its microscopic origin [S]. 

In dilute atomic Bose gas, the quantum theory for the nucleation of Bose- 
Einstein condensation was presented in [6]. This process can be described by 
a single Fokker-Planck equation. Symmetry breaking in the self-evolution of 
strongly non-equilibrium interacting Bose gas can occur only as a result of 
introduction of a small term of the Hamiltonian or as a result of quantum- 
mechanical measurements, which imply non-conservation of particles [7]. 

Dark solitons are one-dimensional phenomena, which in two- or three- 
dimensional (2D, 3D) systems are unstable with a transverse instability. Due 
to the snake-instability [8], they decay into vortex rings [9|. A long- wavelength 
transverse instability that is analogous to matter-waves was observed [10] and 
verified in atomic BECs 1111. 



In harmonic confinement, because of the geometry of the system, a dark 
soliton-like wave moves with different velocities. So, its shape becomes curved, 
which was also observed experimentally [12]. 

An investigation of fractionalized, half-quantum vortices in BEC was 
presented in [13j. They use sodium atoms and rotating optical traps to study 
the creation of individual half-quantum vortices and vortex lattices with the 
help of additional pulsed magnetic trapping potentials. 

They numerically solve the time-dependent coupled GPE of spin-1 and 
use the same damping rate, 7 = 0.03 as we do. Due to the different nature of 
our GPE, our vortices in our presented second case is different than theirs. 

Collisions between two orthogonal solitons and their dynamics were stud- 
ied with coupled GPE in [J3]. It was demonstrated that, the exact vector- 
soliton solution can be obtained with arbitrary tunable interactions. Thus, for 
a sufficiently repulsive interaction in two-component attractive BEC, dark- 
dark solitons can be formed. 

Analytical and numerical investigations of dynamics and modulation of 
ring dark soliton in 2D BEC, were presented in [15] with tunable interaction. 



It was shown that, the stability of a shallow ring dark soliton is stable, when 
the ring is slightly distorted. On the other hand, for large deformation of the 
ring, vortex pairs appear with different dynamical behaviour. 

In this paper, a phase imprinting method [16], [TT] is used to study the for- 
mation of dark solitary waves in two cases, and the formation of the concurrent 
acoustic emissions. 

2 Model 

To explore the instability of soliton-like waves and the density and phase 
profiles of a rotating condensate, we solve numerically the following 2D GPE, 
which governs the time evolution of the order parameter, ^(r, t): 
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Mean field interactions are represented by 52D = V32tt?iu)j— [IE], where u> 



is the angular velocity along the z — axis, I is the magnetic length, l z is the 
characteristic length of the z — axis oscillator and 7 = 0.03 is the dissipation 
modelling the interaction of the system with the surrounding thermal cloud 
[19], [20]. This phenomenological dissipation parameter, 7, was studied in de- 
tail by [21]. Its microscopic justification was introduced by Gardiner et al [22] 
and Penckwitt et al [23]. The chemical potential, /i, can be initially estimated 
in dimensional form as: 



Ngm 

^ = ^y^^> ( 2 ) 

where, N is the number of atoms. The trapping potential becomes: 

V tr (x,y) = ^mu;l(x 2 + y 2 ), (3) 

where, oj± is the radial trap frequency. The quantity is the angular frequency 
of rotation about the z-axis and L z = ih{xd y — yd x ) is the angular momentum 
operator. The density and phase of the condensate together with the superfluid 
local velocity flow are given by p(r,t) = \ip(r,t)\ 2 , S = tan -1 ^^^'^j , and 
v = h/mS7S(r,t). The total energy, E tot is defined as: 

Etot = ^kin + E- mt + E q + Strap, (4) 



where the kinetic energy, -Ekiru the internal energy, E[ nt , the quantum energy, 
E q , and the trap energy, Etrap are given respectively by: 

#kin(t) = / ^ (v/p(x,t)v(x, t)) 2 d 2 r, (5) 

E- mt (t) = Jg(p(^t)) 2 d 2 r, (6) 

E q (t) = J ^ (V ^^j) 2 d\ (7) 

E trap (t) = [ p(x,t)V tT (x)d 2 r. (8) 

The kinetic energy, £7fcj n is a sum of the sound energy, E? soun d , and the vortex 
energy, EVtex UM- 

-C'kin — -C'sound ~r -C'vortcx' yd) 

At a given time, t, the vortex energy, £y or t ex , is achieved by relaxing the 
GPE in imaginary time, which gives the lowest energy state for the vortex 
configuration [23]. Finally, the sound energy is obtained as: -Egound = E^in — 
E'vortox- Dimensionless units are also used, with the units of length, time 
and energy given by: y l /H/(2muj±), u>J and fko±, (harmonic oscillator units, 
h.o.u.). The dimensionless form of the GPE is: 



c \ 9 ^ 



~V 2 + V trap + C\iP\ 2 -p-QL z 



<P, (10) 



where, 



Vfrap = \ {x 2 + V 2 ) , (11) 

and the dimensionless form of the effective coupling constant becomes [261 [271 
EHIEH]: 

and after some transformations, we obtain: 

C = 2V2^^. (13) 

Here, N is the number of atoms, uj_ is the radial trap frequency and l z is the 
characteristic length of the z-axis oscillator. If, C = 1400 and u> 2 S> Clo^/V^tt 



=4> ujI 3> 280ar^. Say, for example uJl > 2800w^ => u z > 53 uj_, where wj_ is 
the radial trap frequency. Throughout this paper, we use C = 1400. 

3 Results 

Dark solitary waves are generated with the help of a phase imprinting 
method, which means a complete absence of fluid (zero density). For the 
wave function an additional phase factor of e ^ is imposed with, Acj) = -k 
(a longitudinal phase change of %). Two cases (A and B) are distinguished 
depending on the presence or absence of the angular frequency of rotation 
about the z — axis, (SI = or £1 ^ 0). 

3.1 Phase imprinting in upper two quadrants for Q = 0. 

We consider a condensate without vortices (0, = 0, no rotation), as 
shown in Fig. [T] (left). The numerical calculations are performed with the 
semi-implicit Crank-Nicholson method in a square box of size between +/— 7 
and +/— 15 (h. o. u.). The origin of the coordinate, (x,y) = (0,0), is at 
the center of the box. The number of grid points vary between 80 and 150. 
Over the duration of the run due to the finite discretization of the numerical 
scheme, the relative change is for example AEt t/Et t — 0.012. 

The phase imprinting condition is: 

ip — ► ip for y > 0, (14) 

and 

ij, -> e i7r V = -^ for y < 0, (15) 

as Fig. H] shows. At t = 200, 7 is set to zero. Fig. [T] (right) presents the 
condensate just after the phase imprinting with large sound waves propagating 
in opposite directions. At t = 200.45 these sound waves reach almost the edge 
of the condensate as we can see in Fig. [2] (left). 

In a confined BEC the density is inhomogeneous. The solitary wave 
embedded in a 2D geometry leads to dominant decay mechanism. Its pertur- 
bation has a longitudinal dynamical instability, and it begins to be excited 
transversely. The applied phase slip of tt accelerates the soliton-like perturba- 
tion due to an interaction with the harmonic trap. 

The slightly curved line becomes torn apart and is accompanied by emis- 
sion of sound waves, which in atomic condensates are density waves. These 
sound waves propagate in opposite direction to the movement of the solitary 



wave and are observed in Fig. [2] (right) and in Fig. [3] (left). Close to the cen- 
ter, the solitary wave propagates faster and close to the boundary slower, see 
Fig. (right) and Fig. [3] (left). 

A system with solitons tend to decay into more stable, lower energy 
structures, namely vortex - anti vortex structures in 2D or vortex rings in 3D. 
Thus, in our case, the soliton-like perturbation bends more and decays into 
pairs of vortices with opposite charges [30j via emission of sound waves, see 
Fig. [3] (left). 

3.2 Phase imprinting in upper two quadrants for Q = 0.8. 

In the second case, the system rotates at angular velocity fi around the 
z-axis. Due to a dynamical instability of the condensate, vortex nucleations 
occur. A stable lattice of 20 vortices is created in a rotating frame. At t = 200, 
the previously presented phase imprinting method is suddenly performed and 
7 is set to zero. The phase imprinting method has an effect on the other 

vortices as well. Initially, there are twenty vortices. After this method, four 
vortex - anti vortex pairs are formed on the line y = 0. Fig. [5] shows the 
displacement of every vortices in the vortex lattice, and the appearance of 
vortex - anti vortex pairs in the phase slip due to the soliton like perturbation. 
The corresponding time for these vortex-motions is: 200 < t < 200.25. 

4 Discussion 

When a discontinuity in the phase is generated, the system tries to 
smooth out this change and generate dark solitary and sound waves. Thus, 
after the application of the suggested method solitary wave is formed in a sys- 
tem with a vortex lattice. This method can be applied to create systems with 
positive and negative vortices in a disc-shaped condensate, whose interactions 
lead to a disordered motion of vortices. 

Fig. [2] (right) and Fig. [3] (left) reveal that, large density oscillations (sound 
waves) are present in the condensate. Due to the confined and finite size of 
the condensate, the sound waves are reflected from the edge of the condensate 
and interact again with the vortex pairs. We can see the anti-correlation 
between E soun d and E VOTtex . When -E S ound is approaching its maximum, E VOTtex 
is approaching its minimum, as shown in Fig. [6] (left) and Fig. [7J We confirm 
that -Esound = -Ekm — -Evortex! compare Fig. [6] (left) and (right) with Fig. [7J 

After applying the phase imprinting method the total energy begins to 
oscillate and later it relaxes to its final value. We remark that, the sound 



energy has the biggest contribution to the change of the total energy, which 
comes from the phase change, from the movement of the vortices and from 
the interaction between vortices and anti - vortices. This is illustrated by the 
Fig. where the phase imprinting method is applied on both y = and x = 
axis (left) or on y = axis (right). Here, E tot is lowered down by 18 units 
to be possible the comparison. From Fig. [9] (left) and (right), we can see the 
real value (in h.o.u.) of the total energy. The compressible kinetic energy is 
directly concerned with the acoustic emission. 

For f2 = 0.8, after a certain transition period, the system enters the state 
that is identical to the original regular vortex lattice. This is explained by the 
fact that although we start with different vortex configurations (no vortices 
versus some artificial vortex - anti vortex pairs), the calculations produce 
identical results. In that case, the negative vortices will move to the boundary 
of the condensate. The corresponding variation of different energies and L z 
are presented by Fig. [U] (right). 

To create turbulence with the double use of this method [3T| on both 
y = and x = axis, immediately after the application of this method, Q is 
set to at t = 200, Fig. [9] (left) shows the energies and the z-component of 
the angular momentum for that case. 

We can compare this method to produce turbulence with a different one, 
when \& = *&* on the upper half part of the box {y > 0). Here, *$>* is the 
complex conjugate of ^(r,t) |29| . In that case, the number of anti - vortices 
and the shape of the condensate are different due to the fact that more anti - 
vortices and annihilations appear. 
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Figure 1: Density profile of the condensate for = 0. Left) at t = 199, before 
the phase imprinting, with no vortices. Right) at t = 200.07, after the phase 
imprinting method with soliton-like wave accompanied by large sound waves 
propagating both in positive and negative directions of the y-axis. 
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Figure 2: Density profile of the condensate corresponding to Fig. [H Left) at 
t = 200.45, the sound waves propagate both upward and downward. Right) 
at t = 200.6, a soliton-like perturbation bends and moves, accompanied by 
sound waves. 
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Figure 3: Density profile of the condensate corresponding to Fig. [TJ Left) at 
t = 200.75, the soliton-like perturbation bends more and decays into vortex 
- anti vortex pairs accompanied by sound waves. Right) at t = 201.15, the 
system decays into vortex pairs. The surviving two pairs from the original five 
pairs of vortices. 
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Figure 4: Phase profile of the condensate at t = 200. Left) two-dimensional 
presentation and Right) three-dimensional presentation. In the random plane 
region, we have large fluctuations of the phase, (the external part of the con- 
densate), due to ip —> 0, Im -ip — > 0, Re ifi — > 0, but the Phase = S = 
tan _1 (Im(V'(r ) t))/Re(V'(r ) t))) ^> 0. 
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Figure 5: 
move in 
vortex - 
200.25). 



For f2 = 0.8, after applying the phase imprinting method the vortices 
a disordered fashion and the soliton-like perturbation decays into 
anti vortex pairs in the center of the condensate, (for 200 < t < 
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Figure 6: For Q = 0.8, Left) Sound energy and Right) Kinetic energy as 
a function of time, from that time when the phase imprinting method was 
applied on y = 0, (here t = corresponds to t = 200). 
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Figure 7: Vortex energy corresponding to FIG. 6. 
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Figure 8: For Q = 0.8, the Total energy —18, Etot — 18, is compared with 
the Sound energy, E soun( i, vs time for Left) the phase imprinting in the upper 
left quadrant and bottom right quadrant (on both y = and x = lines), 
where the Sound energy is the lower /green curve and for Right) the phase 
imprinting in the upper two quadrants (on y = 0), where the Sound energy is 
the lower/red curve. 
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Figure 9: Different energies and L z vs time after the application of the phase 
imprinting method at t = 200 for Left) O, = (the anti - vortices are kept 
inside of the condensate and turbulence is created) and for Right) S7 = 0.8 
(when the system will relax to the initial vortex lattice). Compare L z and 
Eki n in these two cases. 
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